Draft version October 13, 2011 

Preprint typeset using 1^1^^ style eniulatcapj v. 11/26/03 



O 

-)— » 
o 

O 



O 
U 

6 

(N 

> 

m 

o 

l> 

o 



% 



GALAXY STRUCTURE AND MODE OF STAR FORMATION IN THE SFR-MASS PLANE FROM Z - 2.5 TO 

Z-0.1 

Stun Wuyts\ Natascha M. Forster Schreiber\ Arjen van der Wel^, Benjamin Magnelli\ Yicheng Guo^, 

Reinhard Genzel\ Dieter Lutz\ Herve Aussel^, Guillermo Barro'', Stefano Berta\ Antonio Cava'', Javier 

Gracia-Carpio\ Nimish p. Hathi'', Kuang-Han Huang*, Dale D. Kocevski'', Anton M. Koekemoer^, Kyoung-Soo 

Lee"', Emeric Le Floc'h'', Elizabeth J. McGrath'', Raanan Nordon\ Paola PopessoV Francesca Pozzi", Laurie 

Riguccini*, Giulia Rodighiero^^, Amelie Saintonge\ Linda Tacconi^ 

Draft version October 13, 2011 

ABSTRACT 

We analyze the dependence of galaxy structure (size and Sersic index) and mode of star formation 
illsFR and SFRir/SFRuv) on the position of galaxies in the SFR versus Mass diagram. Our 
sample comprises roughly 640000 galaxies at z ^ 0.1, 130000 galaxies at z ^ 1, and 36000 galaxies 
at z ~ 2. Structural measurements for all but the z ^ 0.1 galaxies are based on HST imaging, and 
SFRs are derived using a Herschel-calibrated ladder of SFR indicators. We find that a correlation 
between the structure and stellar population of galaxies (i.e., a 'Hubble sequence') is already in place 
since at least z ^ 2.5. At all epochs, typical star-forming galaxies on the main sequence are well 
approximated by exponential disks, while the profiles of quiescent galaxies are better described by 
de Vaucouleurs profiles. In the upper envelope of the main sequence, the relation between the SFR 
and Sersic index reverses, suggesting a rapid build-up of the central mass concentration in these 
starbursting outliers. We observe quiescent, moderately and highly star- forming systems to co-exist 
over an order of magnitude or more in stellar mass. At each mass and redshift, galaxies on the main 
sequence have the largest size. The rate of size growth correlates with specific SFR, and so does Tjsfr 
at each redshift. A simple model using an empirically determined SF law and metallicity scaling, in 
combination with an assumed geometry for dust and stars is able to relate the observed T^sfr and 
SFRjfi/ SFRuv, provided a more patchy dust geometry is assumed for high- redshift galaxies. 
Subject headings: galaxies: high-redshift - galaxies: stellar content - galaxies: structure 



1. INTRODUCTION 

Deep multi-wavelength lookback surveys carried out 
over the past five years have improved our understand- 
ing of galaxy evolution in the young universe tremen- 
dously. At least over the past 10 Gyr of lookback time, 
a correlation between the rate of star formation and the 
amount of assembled stellar mass in star-forming galax- 
ies (SFGs), dubbed the 'main sequence of star formation' 
(MS), has been observed (Noeske et al. 2007; Elbaz et 
al. 2007; Daddi et al. 2007). While the scatter around 

^ Max-Planck-Institut fiir extraterrestrische Physik, Giessen- 
bachstrasse, D-85748 Garching, Germany 

^ Max-Planck-Institut fiir Astronomie, Konigstuhl 17, D-69117 
Heidelberg, Germany 

^ Astronomy Department, University of Massachusetts, 710 N. 
Pleasant Street, Amherst MA01003, USA 

4 Laboratoire AIM, CEA/DSM-CNRS-Universitc Paris Diderot, 
IRFU/Service d'Astrophysique, Bat. 709, CEA-Saclay, 91191 Gif- 
sur-Yvette Cedex, France 

"■ UCO/Lick Observatory, Department of Astronomy and Astro- 
physics, TJniversity of California, Santa Cruz CA 95064, USA 

° Departamento de Astrofisica, Facultad de CC. Fi'sicas, Univer- 
sidad Complutense de Madrid, E-28040 Madrid, Spain 

"^ Observatories of the Carnegie Institution of Washington, 
Pasadena CA 91101, USA 

* Johns Hopkins University, 3400 N. Charles Street, Baltimore 
MD 21218, USA 

^ Space Telescope Science Institute, 3700 San Martin Drive, Bal- 
timore MD 21218, USA 

^^ Yale Center for Astronomy and Astrophysics, Department of 
Physics, Yale University, New Haven CT 06520, USA 

Dipartimento di Astronomia, Universita di Bologna, via Ran- 
zani 1, 40127 Bologna, Italy 

^^ Dipartimento di Astronomia, Universita di Padova, Vicolo 
deirOsservatorio 3, 35122 Padova, Italy 



the MS is not observed to evolve strongly with redshift, 
the zeropoint does, in the sense that high-redshift SFGs 
form stars at a higher rate than similar mass galaxies 
today. 

Complementing the observational results on the locus 
of SFGs in SFR-Mass space, direct measurements of the 
CO molecular line emission from galaxies residing on the 
MS have revealed molecular gas mass fractions as high 
as 35 (45) % at z ~ 1 (2) (Tacconi et al. 2010). Link- 
ing surface densities of gas and star formation, Genzel et 
al. (2010) recently argued that high-redshift SFGs fol- 
low the same star formation law as local SFGs, unlike 
their merging counterparts that have shorter depletion 
timescales and follow a Kennicutt-Schmidt relation with 
a higher zeropoint (see also Daddi et al. 2010). 

From the observed zeropoint evolution of the MS, 
and tightness of the relation, a picture has emerged in 
which the bulk of SFGs are forming stars gradually over 
timescales that are long (one to two orders of magnitude 
larger) relative to their dynamical times (Genzel et al. 
2010; Wuyts et al. 2011). This is only possible by a con- 
tinuous replenishment of their gas reservoirs. Cold, fila- 
mentary streams of gas have been proposed on the basis 
of cosmological simulations as an efficient mechanism to 
penetrate the galaxies' surrounding halos, and deposit 
new (or recycled) fuel for star formation (Keres et al. 
2005, 2009; Dekel & Birnboim 2006; Dekel et al. 2009). 
Evidence from kinematics strongly supports such a con- 
tinuous (i.e., non- major merger) triggering and mainte- 
nance of star formation for the bulk of SFGs (Genzel et 
al. 2008; Forster Schreiber et al. 2009; Shapiro et al. 



2008; Forster Schreiber et al. 2011a). 

In such a scenario, a secular mode of star forma- 
tion dominates the cosmic SFR history, and the scat- 
ter around the MS can, aside from a contribution by 
measurement uncertainties, be attributed to burstiness 
induced by occasional (minor) merger events. While the 
excess star formation due to mergers may only amount 
to 5-10% (Hopkins et al. 2010; Rodighiero et al. 2011), 
they may play a crucial role in shaping the diversity of 
galaxies observed both at low and high redshift. The 
merger paradigm (Sanders et al. 1988; Hopkins et al. 
2006) links a wide variety of galaxy types in an evolution- 
ary sequence in which the collision between SFGs creates 
tidal torques that channel gas to the center where it trig- 
gers a nuclear starburst (e.g., Barnes & Hernquist 1991, 
1996). When the starburst is accompanied by accretion 
onto the central supermassive black hole(s) (SMBHs), 
this produces AGN emission, and naturally accounts for 
the origin of a correlation between SMBH and bulge mass 
(Di Matteo et al. 2005). Finally, a combination of gas ex- 
haustion, supernova and AGN feedback leads to a rapid 
decline in the SFR, leaving a post-starburst, red and dead 
galaxy as remnant (Wuyts et al. 2010; Snyder et al. 
2011). 

Indeed, such classes of galaxies offset from the MS have 
also been identified observationally out to z ~ 2.5 and 
even beyond: towards the high-SFR end with for in- 
stance bright sub-millimeter galaxies (SMGs, Small et 
al. 1997; Barger et al. 1998; Hughes et al. 1998; see 
Blain et al. 2002 for a review), and towards the low- 
SFR end (Daddi et al. 2005; Kriek et al. 2006; van 
Dokkum et al. 2011). Observations of their structural 
properties, albeit often limited to small samples, could 
lend credence to their merger connection. Several of 
the SMGs imaged in CO show two-component morpholo- 
gies and/or highly disturbed velocity fields, implying a 
major merger nature (Tacconi et al. 2008). Massive 
(~ 10^^ -^0) quiescent galaxies at high redshift are in- 
creasingly more compact relative to similar mass galaxies 
today (fe/re, z=o ~ (1 -l-z)~^'^^°'^) (e.g., van der Wei et 
al. 2008; van Dokkum et al. 2008, but see Mancini et al. 
2010 for a contrasting view), in agreement with expecta- 
tions from binary merger scenarios in which the gas con- 
tent of the progenitors increases with redshift (Khochfar 
& Silk 2006; Hopkins et al. 2009; Wuyts et al. 2010). 
Due to the highly dissipational formation process, Wuyts 
et al. (2010) predict a larger degree of rotation in the 
remnants than in lower redshift early-type galaxies. Ob- 
servational evidence for high-redshift compact quiescent 
disks (McGrath et al. 2008; van der Wei et al. 2011) 
could possibly refiect this, or alternatively may signal 
other quenching mechanisms at play. 

With the high-resolution optical (ACS) and near- 
infrared (WFC3) cameras onboard HST, detailed struc- 
tural measurements can now be obtained for ever larger 
samples, bringing statistical significance also for the more 
extreme populations. In particular, the ongoing Cosmic 
Assembly Near-infrared Deep Extragalactic Legacy Sur- 
vey (CANDELS, Grogin et al. 2011; Koekemoer et al. 
2011) is collecting deep rest-frame optical (and rest-UV) 
space-based imaging over unprecedented areas. At the 
same time, the PACS Evolutionary Probe (PEP, Lutz 
et al. 2011) opens a window on the far-infrared emis- 
sion of large samples of SFGs, and allows the calibration 



of other SFR indicators for use in fields that lack such 
(deep) far-IR data (Wuyts et al. 2011). 

In this paper, we exploit high-resolution imaging and 
a cross-calibrated 'ladder of SFR indicators' to study the 
structure of galaxies and their mode of star formation, 
as function of their position in the SFR-Mass diagram. 
We discuss how properties such as size, surface bright- 
ness profile shape, SFR surface density, and IR/UV ra- 
tios vary along and across the MS, and how these trends 
evolve with redshift (from z '^ 2.5 till today). To this 
end, we assembled samples of unprecedented size, with 
639924 galaxies at z ^ 0.1, 132328 galaxies at z ^ 1, and 
35649 galaxies at z ~ 2. The data quality and sample 
size allow us to extend previous studies of the relation 
between structure and stellar populations in the nearby 
(Kauffmann et al. 2003; Schiminovich et al. 2007; Bell 
2008) and intermediate-redshift (Scarlata et al. 2007; 
Maier et al. 2009) universe, and to build on the ground- 
based explorations of these early epochs (Franx et al. 
2008; Williams et al. 2010). 

The paper is structured as follows. We present an 
overview of the fields and data used, the methods to com- 
pute derived products, and the final sample in Section 2. 
Readers who wish to skip directly to our results, should 
turn to Sections 3 and 4, where we analyze galaxy struc- 
ture and the mode of star formation in the SFR-Mass di- 
agram respectively. Section 5 discusses the observational 
results in a physical context. Finally, we summarize our 
conclusions in Section 6. 

Throughout this paper, we quote magnitudes in the 
AB system, assume a Chabrier (2003) inital mass func- 
tion (IMF), and adopt the following cosmological param- 
eters: (f}M,^A,^) = (0.3,0.7,0.7). 

2. OBSERVATIONS, METHODS AND SAMPLE 

In our analysis, we divide the galaxies in three redshift 
bins: 0.02 < z < 0.2, 0.5 < z < 1.5, and 1.5 < z < 2.5, 
which we will in short refer toasz~0.1,z^l, and 
z ~ 2. The z ~ 0.1 sample is extracted from the Sloan 
Digital Sky Survey (SDSS), as detailed in Section 2.1.1. 
The higher redshift samples are extracted from four deep 
fields that host some of the richest multi- wavelength data 
sets available to date: COSMOS, UDS, GOODS-South, 
and GOODS-North. We describe the data available for 
them in Sections 2.1.2 to 2.1.5, and summarize the meth- 
ods applied to compute derived properties such as pho- 
tometric redshifts, stellar masses, SFRs and structural 
parameters in Section 2.2. 

2.1. Fields and Data 

2.1.1. SDSS + GALEX 

SDSS offers a robust local universe anchor for our study 
of galaxy properties in the SFR-Mass diagram, and its 
evolution with redshift. We use the spectroscopic SDSS 
sample (DR7, Abazajian et al. 2009) in overlap between 
the MPA-JHU and NY-VAGC value-added galaxy cata- 
logs. SFRs and stellar masses were taken from the MPA- 
JHU compilation^^. Briefly, the stellar masses were es- 
timated by fitting Bruzual & Chariot (2003, hereafter 
BC03) stellar population synthesis models with a wide 
range of star formation histories (SFHs) to the broad- 

^'' http;//www. mpa-garching.inpg.de/SDSS/DR7 



band photometry (Salim et al. 2007). Unless explic- 
itly stated otherwise, we use total, dust- and aperture- 
corrected SFRs. SFRs within the SDSS spectroscopic 
fiber were primarily based on Ha, using the calibration 
of Brinchmann et al. (2004), which includes an extinc- 
tion correction based on the Balmer decrement. Aper- 
ture corrections to total use fits to the photometry of the 
outer regions of the galaxies, as described by Salim et 
al. (2007). We use the SDSS (DR7) - GALEX (GR5) 
matched catalog from Bianchi et al. (2011) to com- 
pute the unobscured part of the star formation SFRjjv- 
Sizes and profile shapes were measured by the NY-VAGC 
team^** by means of Sersic fits to the galaxy surface 
brightness profiles (Blanton ct al. 2003, 2005). By de- 
fault, we adopt the structural parameters measured in 
the g band, but we will discuss the dependence on wave- 
length where appropriate. 

2.1.2. COSMOS 

We make use of public multi-wavelength photometry 
from Ilbert et al. (2009) and Gabasch et al. (2008) 
over an effective area of 1.48 deg^ in the COSMOS HST 
field (Scoville et al. 2007; Koekemoer et al. 2007). Our 
working catalog is cut at « < 25 to guarantee sufficient 
signal-to-noise to allow for a reliable determination of 
photometric redshifts and other derived properties. A 
total of 36 medium and broad bands sample the SED 
from GALEX to IRAC wavelengths. At infrared wave- 
lengths, PACS imaging from PEP, extracted using 24 /im 
sources as prior, reaches depths of 5.7 and 11.9 mJy (3cr) 
at 100 /xm and 160 /xm respectively. MIPS 24 /xm imag- 
ing in COSMOS (Sanders et al. 2007; Le Floc'h et al. 
2009) is used down to 60 /iJy. An X-ray catalog based 
on XMM-Newton observations is released by Cappelluti 
et al. (2009). High-resolution J8i4-band imaging with 
HST/ ACS (/814 < 27, 5cr) spanning the entire field, and 
probing the rest-frame B and 2700il. band at z ~ 1 and 
z ~ 2 respectively, forms the foundation for structural 
measurements in COSMOS. 

2.1.3. UDS 

A deep WFC3 -ffigo-selected catalog (-ffieo < 26.7, 5cr) 
with consistent photometry in 16 bands from i? to 8 fiuY 
was constructed on the basis of CANDELS, SEDS (PI 
G. Fazio), SpUDS (PI J. Dunlop), and ancillary ground- 
based data in the ~ 200 arcmin^ CANDELS-covered 
UDS area (Guo et al. in prep). MIPS 24 /im imaging 
to a depth of 30 fiiy (3ct) was carried out as part of the 
SpUDS survey. Ueda et al. (2008) presented the X-ray 
source catalog of the Subaru/XMM- Newton deep survey, 
in which the CANDELS/UDS area is embedded. We per- 
formed the structural measurements on the WFC3 ifieo 
mosaic, which with 4/3 orbits per pointing forms part of 
CANDELS- Wide. The iJigo band probes the rest-frame 
/ and V band at z '-^ 1 and z ~ 2 respectively. 

2.1.4. GOODS-South 

Following similar procedures as for UDS, Guo et al. (in 
prep) constructed a deep WFC3 -ffigo-selected catalog 
(ifieo < 27, 5ct) for the Early Release Science (ERS, PI 
O'Connell) and CANDELS-Deep area in GOODS-South. 

^"^ http://sdss.physics.nyu.edu/vagc 



The catalog contains consistent photometry in 14 pass- 
bands from ?7 to 8 /xm. Spanning the same wavelength 
range, we use FIREWORKS [K^ < 24.3, 5a, Wuyts et 
al. 2008) as supporting multi-wavelength catalog to ex- 
ploit CANDELS- Wide data in the bottom 20% of the 148 
sq. arcmin GOODS-South field. As part of the PEP sur- 
vey (Lutz et al. 2011), PACS photometry was obtained 
to a 5(7 depth of 1.8 mJy, 1.9 mJy and 3.3 mJy at 70 
/xm, 100 /im and 160 /xm respectively, using the position 
of 24 /im sources as prior. The 24 /im imaging itself (30 
/xJy, 5ct) was taken from Magnelli et al. (2009). Where 
appropriate, we select X-ray sources from the Chandra 2 
Ms source catalog by Luo et al. (2008). We performed 
the structural measurements on WFC3 ifieo-band imag- 
ing drizzled to a fffOG/pixel scale, obtained as part of 
the ERS and CANDELS programs. The WFC3/IR ob- 
servations of the ERS were obtained in F098M, F125W 
and F160W for a total depth of 2 orbits per pointing 
in each filter (Windhorst et al. 2011) and the mosaics 
used here were produced using MultiDrizzle (Koekemoer 
et al. 2002) following the approach outlined in Koeke- 
moer et al. (2011). The 4-epoch CANDELS-Deep area 
is similar in depth, while CANDELS- Wide received half 
the integration time. 

2.1.5. GOODS-North 

We derived photometric redshifts and stellar masses for 
galaxies in GOODS-North based on a z -I- X-selected cat- 
alog with photometry in 16 bands from GALEX to IRAC 
wavelengths (Berta et al. 2010). Additional longer wave- 
length information at 24 /im (30 /iJy, 5cr), 100 /im (5.1 
mJy, 5(t), and 160 /xm (8.7 mJy, 5a) was acquired as part 
of the GOODS and PEP surveys. Where appropriate, we 
select X-ray sources from the Chandra 2 Ms source cata- 
log by Alexander et al. (2003). Structural measurements 
in GOODS-North are based on the GOODS ACS zgso- 
band imaging (Giavalisco et al. 2004), corresponding to 
the rest-frame B and 3000A band at z ^ 1 and z ^ 2 
respectively. . 

2.2. Derived Galaxy Properties 

For all deep fields, we use identical procedures to derive 
photometric redshifts, stellar masses, SFRs, and struc- 
tural parameters, hence optimizing the consistency of the 
combined data set. 

2.2.1. Photometric Redshifts 

We determined photometric redshifts by fitting a su- 
perposition of 6 galaxy templates to the observed broad- 
band SEDs using the photometric redshift code EAZY 
(Brammer et al. 2008). For the 6955 out of 132328 
(5.3%) of galaxies at 0.5 < z < 1.5, and 497 out of 35649 
(1.4%) of galaxies at 1.5 < z < 2.5 that have a spec- 
troscopic redshift from one of the many spectroscopic 
campaigns in the above four fields (see, e.g., Vanzclla et 
al. 2008; Lilly ct al. 2009; Barger et al. 2008), we adopt 
the Zspec as redshift in our analysis. 

Comparing the photometric redshifts Zphot and spec- 
troscopic redshifts Zgpec of spectroscopically confirmed 
galaxies, we find a median and scatter (normalized me- 
dian absolute deviation) in Az/(1 + z) of (-0.001; 0.015) 
in the z ~ 1 and (-0.007; 0.052) in the z ~ 2 red- 
shift bin. Naturally, these uncertainties only apply to 



galaxies with similarly bright magnitudes as represented 
in the spectroscopic sample based on which the scat- 
ter was computed. E.g., the public zCOSMOS-bright 
spectroscopic survey (Lilly et al. 2009) reaches down to 
i = 22.5. For this spectroscopic sample in COSMOS, we 
measure a scatter in Az/(1 + z) of 0.010. As detailed 
in Appendix B.3, the estimated uncertainty increases to 
aNMADiAz/{l + z)) ^ 0.048 for COSMOS galaxies with 
24 < z < 25, which account for 40% of the galaxies in 
that field, and 30% of our entire sample. The estimated 
uncertainty is consistent with earlier work by Ilbert et 
al. (2009). For the other deep fields as weU, the increase 
in the photometric redshift uncertainty with respect to 
the spectroscopic sample amounts to a factor of ~ 5 for 
the faintest galaxies that enter our analysis. An exten- 
sive analysis in Appendix B.3 shows that our results are 
robust against photometric redshift uncertainties, even 
when accounting for the fact that the Zphot quality of the 
entire sample is likely lower than that of the (typically 
brighter) spectroscopically confirmed subsample. 

2.2.2. Stellar Masses 

We used FAST (Kriek et al. 2009a) to fit BC03 mod- 
els with exponentially declining SFHs to the Xobs < 8 /^m 
broad-band SEDs of the galaxies. Following the recipe 
presented by Wuyts et al. (2011), we impose a min- 
imum e-folding time of 300 Myr. The time since the 
onset of star formation was allowed to vary between 50 
Myr and the age of the universe at the observed redshift. 
Our choice of SFHs guarantees consistency with most of 
the literature, and in addition has the advantage (as op- 
posed to, e.g., only exponentially increasing SFHs) that 
it is also applicable to lower redshift galaxies and high- 
redshift quiescent galaxies, that clearly formed stars at 
a higher rate in their past than at the epoch of obser- 
vation. Maraston et al. (2010) showed that adopting 
exponentially increasing SFHs for star-forming galaxies 
at z ~ 2 leads to negligible changes in the estimated 
stellar masses, unless also a stringent constraint on for- 
mation redshift is imposed, in which case the estimated 
stellar mass of these galaxies can increase by a few 0.1 
dex. Solar metallicities were assumed throughout, and 
dust attenuation was modeled using the Calzetti et al. 
(2000) reddening law, with visual extinctions in the range 
< Ay < 4. For a detailed description of the depen- 
dence of estimated stellar masses on the assumed metal- 
licity and attenuation law, we refer the reader to Wuyts 
et al. (2007). We revisit the assumption of solar metal- 
licity in Section 5.2 in the context of a model for IR/UV 
ratios that involves a metallicity scaling, and conclude 
that it has a negligible impact on our results. As for the 
SDSS galaxies, we assumed a universal Chabrier (2003) 
IMF for all galaxies in the lookback surveys. 

2.2.3. Star Formation Rates 

Exploiting deep PACS photometry in the GOODS- 
South field, Wuyts et al. (2011) established a conti- 
nuity across SFR indicators, ranging from UV -I- FIR, 
over UV + MIR, to SED modeled SFRs. We adopt this 
cross-calibrated 'ladder of SFR indicators' to obtain an 
estimate of the SFR for each galaxy in our sample indi- 
vidually. Briefiy, if the galaxy is detected in an IR band, 
we compute the total SFR by summing the unobscurcd 



and re-emitted emission from young stars, following Ken- 
nicutt (1998): 

SFRuv+iR [Moyr-i] = 1.09x10-^° (L/K+3.3i280o)/io 

where L2800 = i'Li,{2800A) was computed with EAZY 
from the best-fitting SED. Here, the total IR luminos- 
ity Ljji = L(8 — lOOO/irn) is derived monochromatically 
from the longest wavelength IR band that is significantly 
(> 3ct) detected: either PACS 160, 100, or 70 fim, or 
MIPS 24 ^m. We use the mid- to far-infrared SED tem- 
plate by Wuyts et al. (2008) for the conversion to Lju. 
This conversion leads to a consistency between 24 fim 
and PACS-derived L/ij, unlike locally calibrated tem- 
plate sets that have been proven to overestimate the Ljji 
of z > 2 galaxies based on their rest-frame 8 /im emis- 
sion (Nordon et al. 2010; Elbaz et al. 2010; Wuyts et al. 
2011). 

For galaxies that lack an IR detection, we adopt the 
SFR associated to the best-fit stellar population synthe- 
sis model from Section 2.2.2. For the low- to intermediate 
SFR regime, where SFRjjv+iR is only attainable in the 
deepest of our four fields (GOODS-South), Wuyts et al. 
(2011) demonstrated that SED modeled SFRs are con- 
sistent with those obtained from Equation 1. As for the 
stellar masses, the values of SFRsed, obtained without 
imposing explicit constraints on the formation redshift, 
would not change significantly when assuming exponen- 
tially increasing rather than decreasing SFHs (Maraston 
et al. 2010). 

2.2.4. Structural Parameters 

In this paper, we limit our analysis of galaxy structure 
to a parametric approach, focussing on two basic prop- 
erties of the surface brightness distribution: its extent 
and cuspiness. The former is expressed as the circular- 
ized effective radius rp, = ?'p..m.o.ior x \/b/a correspond- 
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ing to the best-fit Sersic profile. The latter is quanti- 
fied by the Sersic index n, with n restricted to the range 
0.2 < n < 8. We adopt the structural measurements per- 
formed on the longest wavelength high-resolution imag- 
ing available for each field (WFC3 for GOODS-S and 
UDS, ACS for GOODS-N and COSMOS, see Section 
2.1). The GALAPAGOS package (Harden et al. 2009) 
or equivalent procedures were used to pre-detcrminc the 
sky background and the neighboring sources that need to 
be simultaneously fit or masked, and to automate fitting 
Sersic profiles to the 2D surface brightness distributions 
with GALFIT (C. Y. Peng et al. 2010). GALFIT takes 
into account convolution by the point spread function 
(PSF). 

Since the longest wavelength at which high-resolution 
imaging is presently available varies between our four 
deep fields, our study lacks the potential to consistently 
probe the rest-frame optical regime, and fully control 
morphological k-corrections, if present. However, since 
we are interested not only in galaxies on the main se- 
quence, but also in those rare systems that form the 
high-SFR tail, we prefer to sacrifice this aspect of our 
analysis, rather than area. As a sanity check, we com- 
pared sizes and Sersic indices measured in the zgso and 
Hieo imaging of GOODS-S (see Appendix B.2). We find 
the median ACS - WFC3 deviations to be small (-0.01 
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Fig. 1. — Surface brightness profile sfiape in the SFR-Mass diagram. A 'structural main sequence' is clearly present at all observed 
epochs, and well approximated by a constant slope of 1 and a zeropoint that increases with lookback time (white line). While SFGs on 
the MS are well characterized by exponential disks, quiescent galaxies at all epochs are better described by De Vaucouleurs profiles. Those 
galaxies that occupy the tip and upper envelope of the MS also have cuspicr light profiles, intermediate between MS galaxies and red & 
dead systems. 



in AlogTe and -0.10 in n at 1.5 < z < 2.5) relative 
to the trends that will be discussed in this paper, sug- 
gesting that any biases from morphological k-corrections 
have only a limited impact on our conclusions. Bond et 
al. (2011) arrive at a similar conclusion comparing rest- 
frame optical and rest-frame UV size measurements of 
z ^ 2 galaxies. Moreover, we also verified that the re- 
sults obtained for each of the four fields individually are 
consistent with those presented for the combined data set 
in this paper, albeit with more noise due to the smaller 
number statistics (see Appendix B.l). 

While beyond the scope of this paper, the wavelength 
dependence of structural properties can potentially re- 
veal interesting clues on the physical processes shaping 
galaxies in the young universe (see, e.g., Guo et al. 2011; 
Szomoru et al. 2011; and Cassata et al. 2011, who find 
evidence for mild morphological k-corrections such that 
the centers of galaxies in their samples tend to be some- 
what redder). For an in depth discussion of resolved stel- 
lar populations inferred from spatial variations in color, 
we defer the reader to Wuyts et al. (in prep). 

2.3. Sample Selection 

Our four deep fields are not uniform in depth. Con- 
sequently, they have different completeness limits in the 
SFR-Mass diagram. Since we are mainly interested in 
individual galaxy characteristics rather than abundances 
(number or mass densities) as function of position along 
or across the MS, we refrain from applying any incom- 
pleteness corrections to the observed populations. Under 
the premise that, at any given redshift, galaxies form a 
two-parameter family described by their mass and SFR, 
this approach works well. I.e., a median galaxy property 
can be computed reliably, and unbiased by any complete- 
ness issues, based on the objects observed in a given bin 
of SFR-Mass space. 

Our final sample comprises 639924 galaxies at 0.02 < 
z < 0.2, 132328 galaxies at 0.5 < z < 1.5, and 35649 



galaxies at 1.5 < z < 2.5. The relative breakdown in 
galaxies of different masses is determined by the depth 
of the observations, and the stellar mass function at the 
respective redshifts. Above M > 10^° Mq our sample 
counts 532131, 31127, and 8895 galaxies at z '^ 0.1, z ^ 
1, and z ^ 2 respectively. Above M > 10^^ Mq, the 
numbers drop to 147922, 2767, and 1059 galaxies at z ~ 
0.1, z '^ 1, and z ~ 2 respectively. An overview of the 
sample size per field is provided in Table 1. 

3. RESULTS ON GALAXY STRUCTURE 

3.1. Profile Shape 

We start by analyzing the surface brightness profile 
shape as a function of position in the SFR-Mass diagram 
in Figure 1. The three panels show from left to right 
the z ~ 0.1, z ~ 1, and z '^ 2 bins respectively. In- 
stead of indicating the relative abundance of galaxies in 
different regions of the diagram, we use the color-coding 
to mark the median value of the Sersic index n of all 
galaxies in each [SFR, M] bin. For displaying purposes, 
we restrict the range of the colorbar to 1 < n < 4, and 
assign the same color as n = 1 and n = 4 to bins with 
median rt < 1 or median n > 4 respectively. The frac- 
tions (/n<i; fn>i) of galaxies lying outside these bounds 
amounts to (0.09; 0.24), (0.41; 0.14), and (0.41; 0.16) at 
z ^ 0.1, z ^ 1, and z ~ 2 respectively. The fraction 
of [SFR, M] bins with median n outside this range is 
small: (0.02; 0.11) at z - 0.1, (0.14; 0.15) at z ~ 1, and 
(0.11; 0.11) at z ~ 2. The resulting diagrams present a 
remarkably smooth variation in the typical galaxy pro- 
file shape across the diagram. Moreover, despite the loss 
of information on number densities, the so-called main 
sequence of star formation is immediately apparent, and 
its presence persists out to the highest observed redshifts. 
This 'structural main sequence' consists of galaxies with 
near-exponential profiles (n sa 1) and shows a similar 
behavior as the conventional 'number main sequence' as 
identified on the basis of number densities in the SFR- 
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Fig. 2. — Top panels: SFR histograms relative to the 'structural main sequence' for galaxies in the mass range 5 X 10 Mq < M < 10 Mq . 
Bottom panels: Sersic index as function of deviation from the MS, measured along the SFR axis, for the same 10^'^ — 10^^ Mq mass slice. 
The black curve and error bars indicate the median in bins of Alog{S FR), and respective errors in the median. The grey polygon marks 
the central 68th percentile of the distribution. Colored dotted curves in the left-hand panel illustrate how the 2 ~ 0.1 relation changes 
when adopting structural measurements performed at shorter/longer wavelengths than the SDSS g band. Relative to the MS, a correlation 
between profile shape and star-forming activity is present and looks similar at all epochs. The relation is not monotonic, but rather shows 
a reversal at the high-SFR end. 



Mass diagram (e.g., Noeske et al. 2007; Elbaz et al. 
2007; Daddi et al. 2007). Namely, an upward shift of the 
zeropoint is observed with increasing lookback time. At 
each epoch, the MS in Figure 1 is well approximated by a 
slope of unity {white line). The SFR at which the median 
n reaches a minimum in a mass slice around log(M) — 10 
roughly coincides with the mode of the \og{SFR) distri- 
bution in that mass slice, but depending on the fitting 
method and sample definition used to weed out quiescent 
galaxies, a somewhat shallower slope than unity may be 
measured for the 'number main sequence' at the massive 
end (see, e.g., Rodighiero et al. 2010). 

Below the structural MS, a cloud of galaxies with 
cuspy, near-De Vaucouleurs (n ~ 4) profiles is visible. 
This population of massive quiescent galaxies is present 
at all observed epochs. Our first and foremost conclusion 
from Figure 1 is therefore that a correlation between the 
structure and stellar populations of galaxies is already 
in place since z ~ 2.5. In the local universe, the pres- 
ence of such a correlation has been described in depth 
by, e.g., Kauffmann et al. (2003) and Brinchmann et 
al. (2004). Scarlata et al. (2007) demonstrated that 
photometrically selected massive quiescent galaxies are 
already dynamically relaxed (i.e., have the morpholog- 
ical appearance of early- type galaxies) at z ^ 1. Also 
out to z ^ 1, Maier et al. (2009) exploited the spectro- 
scopic zCOSMOS survey to reveal a correlation between 
the specific SFR of galaxies and the Sersic index of their 
surface brightness profiles. At larger lookback times, our 
results are in agreement with previous findings, based on 
smaller samples, reporting a similarly early emergence of 
the Hubble sequence on the basis of stellar surface mass 
density, galaxy size and visual appearance (Franx et al. 
2008; Toft et al. 2009; Kriek et al. 2009b). We find that, 
across cosmic time, the typical Sersic index of galaxies 
is not optimally described as a function of their abso- 



lute SFR, or even their absolute specific SFR, but rather 
as a function of their position relative to the MS at the 
epoch of observation. The correspondence between mass, 
SFR and structure, as quantified by the Sersic index, is 
equivalent to the Hubble sequence. Based on samples 
of unprecedented size, we see that such a sequence al- 
ready existed at z ~ 2; bulge-dominated morphologies 
go hand in hand with a more quiescent nature. First in- 
dications of such a correlation between Sersic index and 
specific SFR out to high redshift were recently reported 
by Szomoru et al. (2011) based on 27 galaxies at z '-^ 1 
and 16 galaxies at z '--^ 2. 

In detail, an additional interesting aspect is revealed 
by the n distribution in SFR-Mass space. The sur- 
face brightness profile shape does not vary monotonically 
across the main sequence, but instead shows a reversal 
in the high-SFR tail of the distribution, such that the 
upper envelope of the MS is composed of galaxies whose 
typical cuspiness is intermediate between those of normal 
MS galaxies and the quiescent population below the MS. 
This is illustrated further in Figure 2, where we plot a cut 
through of the 5 x 10^ Mq < M < 10" Mq mass slice. 
The top panels show the histogram of deviations from the 
MS measured along the SFR axis: Alog(S'Fi?). The top 
left panel illustrates the well-known bimodality between 
star-forming and passive galaxies in the local universe 
(e.g., Kauffmann et al. 2003). A similar bimodality, al- 
beit with a reduced amplitude of the quiescent peak, is 
visible at z ~ 1 (and becomes more pronounced if we 
were to limit ourselves to higher masses only). While 
both populations are also found in our highest redshift 
bin, the quiescent fraction is further reduced. Incom- 
pleteness in our shallowest field, as well as uncertainties 
in photometric redshifts and estimated SFRs inhibit the 
identification of a true bimodality in the z '^ 2 histogram. 

The bottom panels of Figure 2 show the median de- 
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Fig. 3. — Galaxy size in the SFR-Mass diagram. Moving up the main sequence, towards the high-mass end, galaxies grow bigger (i.e., 
a size- mass relation is observed at all epochs). With respect to both quiescent galaxies, and also the high-SFR tail, galaxies on the MS 
have the largest size. Finally, in any given [log(A/); log(5-FiJ)] bin, galaxies grow with time. Quantitatively, this growth is strongest for 
quiescent (low SFR/M) galaxies. 



pendence of n on Alog(5Fi?) {black curve), as well as 
errors on the determination of the median [black error 
bars). With grey polygons, we mark the central 68th 
percentile. Clearly, a variety of profile shapes is mea- 
sured, even at a given Alog(S'Fi?), and particularly in 
the z ~ 1 and z ~ 2 bins. This stems partly from mea- 
surement uncertainties, that dominate in the low sur- 
face brightness regime, but in addition likely means that 
in detail galaxies are more complex than a simple two- 
parameter {[SFR,M]) family. Interestingly, the scatter 
seems to be minimized for galaxies that reside on the 
MS. Forster Schreiber et al. (2011) present a detailed 
analysis of large and clumpy MS galaxies at z ~ 2, that 
exhibit shallow surface brightness profiles (n < 1). Such 
a population is also present in our sample (see Figure 2 
and the above quoted /„<i). Focussing on the overall 
trend, we find the results at z ^ 1 and z ~ 2 to be very 
similar to those established at z ^ 0.1. Namely, galaxies 
on the MS resemble exponential disks, high-SFR outliers 
have a centrally enhanced surface brightness profile, and 
once galaxies end up below the MS, they quickly reach 
a plateau of n w 4. We discuss the physical implications 
of this empirical relation in Section 5.1. 

We verified that, when excluding X-ray sources from 
our analysis, the same median dependence of the sur- 
face brightness profile shape on the location of galaxies 
in the SFR-Mass diagram remains present. This is the 
case even for the fields with the deepest X-ray data. The 
same robustness also applies to the other galaxy proper- 
ties discussed later in this paper. Our test indicates that 
the observed trends (particularly the enhanced Sersic in- 
dices of the high-SFR outliers at a given mass) is not 
merely due to a nuclear, non-stellar point source. When 
focussing on the X-ray detected sources alone, we find 
that they occupy the mass regime above logM > 10, 
and span a wide range in SFRs, from the quiescent class 
to on and above the MS. Within a given bin, their sur- 
face brightness profiles differ in the sense that they have 
higher Sersic indices and smaller sizes than normal galax- 



ies of the same [SFR, M] . Point source contamination 
would drive the structural parameters in the observed di- 
rection, but extensive simulations (see, e.g., Simmons & 
Urry 2008) are required to constrain whether in addition 
there is any evidence for a different stellar structure of 
the hosts. 

Before moving to the next structural property, we have 
a closer look at the z ^ 0.1 bin in Figure 2, where we use 
structural measurements in five passbands to investigate 
how the observed trend depends on wavelength. The 
dotted colored curves in the bottom left panel of Figure 
2 illustrate the median relation derived from Sersic fits in 
the SDSS u, r, i, and z bands. Together with the default 
g-band profile fits (solid black curve), a trend emerges in 
which the overall relation between n and A \og{SFR) is 
remarkably robust against morphological k-corrections. 
The largest difference is observed for galaxies that lie on 
the MS {Alog{S FR) = 0), such that a cuspier profile 
is inferred from redder waveband imaging. This comes 
as no surprise, as it is well known that many normal 
SFGs in the nearby universe are composed of a (red, 
n « 4) bulge and (blue, n w 1) disk component^^. At low 
and high values of A\og{SFR), the Sersic index levels 
of at similar values (n = 4 and n = 2 respectively), 
irrespective of wavelength. This suggests that the light 
concentration in the centers of these galaxies corresponds 
to a mass concentration, and is not due (only) to the large 
luminosity of a young stellar population. 

3.2. Size 

Having established the existence of a 'Hubble sequence' 
out to z ^ 2.5, we now turn to the zeroth order structural 
measurement of size, and address how it varies along and 
across the MS at each of our observed epochs. Figure 3 
is similar in nature to Figure 1, but the color-coding is 
now used to mark differences in the spatial extent of the 

^^ The relative fraction of such bulged disks is observed to de- 
crease with redshift, at least out to ^ ~ 1 (Sargent et al. 2007). 
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Fig. 4. — Surface density of star formation (Egj?^) in the SFR-Mass diagram. Lines of equal Hgpji run approximately parallel to the 
MS, as do lines of constant SFR/M. 



surface brightness distribution (re). As in Figure 1, a 
qualitatively similar behavior is observed at all redshifts. 
Overall, more massive galaxies tend to be larger. How- 
ever, the mass-size relation is not fundamental in the 
sense that at a given mass substantial variations in size 
occur that are not random, but instead correlate with 
the offset from the MS (see also the Appendix, Figure 
8). This translates to different mass-size relations fol- 
lowed by passive and star-forming galaxies, as reported 
previously by Shen et al. (2003) based on SDSS, but 
also by Williams et al. (2010) and Szomoru et al. (2011) 
who extended the analysis for massive galaxies out to 
z = 2, based on ground-based K-band imaging in UDS 
and WFC3 iJieo-band imaging in the Hubble Ultra Deep 
Field respectively. 

Consistent with previous results by, e.g., Franx et al. 
(2008) and Toft et al. (2009), high-redshift quiescent 
galaxies are more compact than the bulk of SFGs at the 
same epoch. Also at z ^ 0.1, galaxies below the MS (at 
least at M < 10^^ Af©) tend to be smaller than those on 
the MS. The high-SFR tail of the z -- 0.1 distribution 
differs from normal MS galaxies in the sense that they 
tend to have higher surface mass densities (i.e., are more 
compact for a given mass). As we march up in redshift, 
our measurements are consistent with a similar behavior 
across all lookback times considered. 

Comparing the panels for the three redshift bins, it 
is immediately apparent that galaxies grow over time. 
We quantified the growth by fitting a rpXz)/rf,{z = 0) = 
(1 4- z)" relation to the galaxies in each [SFR,M] bin. 
It is important to note that this reflects the size evolu- 
tion of the galaxy population as a whole, and not per 
se that of any given galaxy individually, as galaxies will 
move from one [SFR, M] bin to another as they build 
up stellar mass, and new systems will be added over 
time. We find a to be varying across the SFR-Mass di- 
agram, from a = -l.2toi below \og{SFR/M) < -11 
to a = -0.4J:o:7 above \og{SFR/M) > -11. Here, the 
range indicates the normalized median deviation over the 
different [SFR, M] bins. Splitting the latter class into 



galaxies above and below the MS at z ^ 0.1 (white line in 
the left-hand panel of Figure 3), we obtain a = — 0.2_Qg 

and a = — O.Q^pg respectively. In other words, the rate 
of size evolution correlates with the specific star forma- 
tion rate, reaching the lowest values in the (especially 
massive) quiescent population. 

4. RESULTS ON STAR FORMATION MODE 

Having established the variation of galaxy structure 
along and across the MS, we now do the same for two 
parameters that characterize the mode of star formation: 
its surface density (Section 4.1), and its breakdown into 
unobscured and dust-enshrouded star formation (Section 
4.2). A discussion of the physical implications of the 
observed relationships follows in Section 5.2. 

4.1. Surface Density of Star Formation 

Figure 4 presents the dependence of the SFR surface 
density 



EsFR = SFR/2Tirl 



(2) 



on the position of galaxies in the SFR-Mass diagram. 
By construction, this Figure has the same information 
content as Figure 3, but its interpretation in terms of 
mode of star formation is more straightforward. Galaxies 
that form stars more actively are not just upscaled ver- 
sions of lower SFR systems. Their ISM conditions differ 
in the sense that also normalized by area they form more 
stars per unit time. Lines of constant ^sfb. I'un diag- 
onally, implying a tighter correlation with specific SFR 
than with SFR or stellar mass separately. Schiminovich 
et al. (2007) discussed this trend extensively for galax- 
ies in the nearby universe. At higher redshifts, this is 
consistent with recent findings based on rest-frame UV 
size measurements by Elbaz et al. (2011). We confirm 
that the same trend is seen when rest-frame optical size 
measurements are used instead. 

We contrast 'Ssfr to SFR/M in Figure 5, where we 
plot the 25th and 75th percentile contours {solid and 
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Fig. 5. — Surface density of star formation as function of spe- 
cific star formation rate. Solid and daslied contours comprise 25% 
and 75% of the galaxies respectively. At all observed epochs, 
galaxies line up along the same fundamental relation, but towards 
higher redshifts the bulk of galaxies shifts to the higher T,sfr (or 
SFR/M) end of the relation. 



dashed lines) for each of the considered redshift bins. 
Over 10 Gyr of lookback time, the bulk of galaxies fol- 
low a similar, linear relation between specific SFR and 
surface density of star formation. As we probe higher 
redshifts, the upper end of this relation becomes more 
densily populated. This trend can naturally be explained 
by increased gas mass fractions in the young universe 
(e.g., Tacconiet al. 2010). 

In equation 2, we adopted the same size measurement 
as discussed in Section 2.2.4 and 3.2 to represent the ra- 
dius within which half the star formation takes place. 
I.e., we adopted the best- fit r^ of a single Sersic profile 
fit. At z ~ 0.1, we used the g band by default. For the 
higher redshift samples, we used the longest wavelength 
high-resolution image available (i^igo in the case of UDS 
and GOODS-S, zsso and hu for GOODS-N and COS- 
MOS respectively). In principle, radial variations in SFH 
and/or dust attenuation could bias the measurement of 
a half-SFR radius. E.g., in the case of a superposition 
of a star-forming disk and a red and dead bulge, the to- 
tal half-light radius would contain less than half the star 
formation (particularly when measured in the red). In 
contrast, in the case of a galaxy where old and young 
stars trace the same distribution, but the dust column is 
centrally enhanced, the observed r^, would contain more 
than half of the star formation (particularly when mea- 
sured in the blue). We take two simple approaches to 
investigate the impact of these potential biases on our 
results. First, we repeated our analysis at z '--^ 0.1 adopt- 
ing each of the ugriz SDSS bands to determine rg and 
based thereupon Y,sfr- We find that this induces a mi- 
nor change only, in the sense that, for galaxies below the 
main sequence, T,sfr increases by a few 0.1 dex from 
the bluest to the reddest band (see the left-hand panel 
of Figure 8 in the Appendix). Overall, the pattern of 
'^SFR in SFR-Mass space does not depend appreciably 
on the adopted photometric band at z ^ 0.1. Second, we 
stress that at z ~ 1 as well as at z ~ 2, the dependence 
of TisFR on location in the SFR-Mass diagram remains 



the same (although with somewhat enhanced noise due 
to the smaller number statistics) when considering the 
fields with rest-frame optical and rest-frame UV imaging 
separately (see Appendix B.l). 

4.2. Obscuration of Star Formation 

Finally, we consider the ratio of re-emitted to 
unobscured star formation. What we denote as 
'SFRiFtJSFRuv' is for IR-detected galaxies computed 
as the ratio of the IR and UV contribution to the total 
SFRjjv+iR (see Equation 1): L/ij;/3.3i2800- For those 
galaxies in the deep lookback surveys that lack an IR de- 
tection, and for which we adopted the best-fit SFR from 
stellar population modeling, we derive SFRirt/ SFRjjv 
as 10°''^28oo _ 1 Here, A2800 is the extinction at 2800A 
as inferred from the SED modeling. This approach has 
the advantage of always producing positive values, un- 
like adopting [SFRsed — S F Rjjv) / S F Ruv as defini- 
tion, which can reach negative values when SFRjjy > 
SFRsed- The latter situation can occur for quiescent 
galaxies with low ratios of present to past-averaged SFR, 
in which case the rest-frame 2800A emission contains a 
substantial contribution from older stars. 

For the z '^ 0.1 sample, total SFRs (based on dust- and 
aperture-corrected Ha measurements, see Section 2.1.1) 
are available from the MPA-JHU release, and we com- 
puted SFRjjv as for the higher redshift galaxies, where 
we interpolated L2800 from the GALEX NUV and SDSS 
u-band photometry. The derived SFRm/SFRuy = 
[SFRtot - SFRuv)/SFRuv at z - 0.1 are well defined 
for the star-forming population, but ill-constrained for 
the passive population for the above reasons. We there- 
fore choose to exclude the passive z '■^ 0.1 galaxies from 
our analysis. 

The resulting SFR-Mass diagrams, color-coded by 
SFRjji/SFRuv: are presented in Figure 6. Again, the 
qualitative appearance of the diagrams looks similar over 
the wide range of redshifts considered. At a given mass, a 
proportionally larger fraction of the star formation activ- 
ity is revealed by reprocessed IR emission as we consider 
more actively star-forming systems. A positive correla- 
tion between bolometric luminosity and obscuration was 
presented previously for galaxies in the nearby universe 
(Heckman et al. 1998) and at z ~ 2 (Reddy et al. 2010), 
but these authors did not trace the trend with SFR while 
specifically controlling for stellar mass. Likewise, an in- 
crease in SFRin/SFRijv is observed as we move along 
the MS from the low- to high-mass end. Such a trend was 
not seen as clearly in T^sfr (Figure 4), a finding that we 
will interpret as due to metallicity variations across the 
diagram in Section 5.2. Finally, objects forming the up- 
per rim of the galaxy distribution, also at lower masses, 
show enhanced IR emission. The latter trend is most 
prominent in the z ^ 1 bin, but a hint of the same fea- 
ture is seen at both lower and higher redshift as well. 

5. DISCUSSION 

5.1. The Relation between Peak Starbursts and 
Quiescent Galaxies 

The observation that a precursor of the Hubble se- 
quence, with disk-like SFGs and cuspier quiescent sys- 
tems, was already in place 10 Gyr ago (see Section 3.1) 
has important implications. To name one, a mechanism 
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Fig. 6. — Ratio of re-emitted over unobscurcd star formation in the SFR-Mass diagram. The highest SFRjji/ SFRuy ratios are found 
at the upper rim, and at the high-mass end of the MS. 



for galaxies to transition from the star-forming to qui- 
escent phase must be universally present across cosmic 
time, at least out to z ~ 2.5. It is during this epoch 
around z ~ 2 that the cosmic SFR density is observed to 
reach its peak (e.g., Hopkins & Beacom 2006), and that 
large accretion rates onto galaxy halos are expected from 
cosmological hydrodynamical simulations (Keres et al. 
2005, 2009; Dekel & Birnboim 2006; Dekel et al. 2009). 
We can therefore exclude gas exhaustion as sole driver of 
the transition from the active to passive phase. Instead, 
a quenching mechanism must be at play. Moreover, the 
observed relation between structure and stellar popula- 
tions implies a causal connection between the quenching 
and the morphological transition. 

As early as by Toomre & Toomre (1972), mergers 
have been invoked as a potential mechanism to transi- 
tion from late to early Hubble types. In such a scenario, 
much of the gas content of the progenitors is consumed 
by star formation in the nucleus, building up a central 
cusp. Violent relaxation distributes the stars that were 
already formed prior to final coalescence in a near-De 
Vaucouleurs [n — 4) distribution. Feedback from short- 
lived massive stars that explode as supernovae, as well 
as energy that is released from accreting material onto 
the SMBH(s) and that couples to the surrounding ISM 
drives out the remaining gas and heats up the galaxy 
halo preventing new infall (e.g., Hopkins et al. 2006, and 
references therein). 

In this light, the reversal of the decreasing trend of n 
with A \og{SFR) (see Figure 2) may signal that we are 
observing the galaxies at the high-SFR tail of the dis- 
tribution during the peak of their SFH, in the process 
of building up a central concentration of stellar mass, 
much like that typically observed in the quiescent pop- 
ulation. This would also explain their increased SFR 
surface densities and small sizes relative to normal SFGs 
on the MS. While galaxies above and below the MS by 
definition form the opposite extremes in terms of SFR, 
they may thus share a more similar nature in terms of 



structural shape. According to the merger evolutionary 
sequence, this would place them as nuclear starbursts in 
between MS galaxies and quiescent remnants. 

A crude duty cycle argument reveals that, if all qui- 
escent galaxies at 0.5 < z < 1.5 with 10^° Mq < M < 
10^^ A'Iq and Alog(iS'Fi?) < — 1 underwent a starburst 
of A log(S'Fi?) > 0.5 prior to quenching, the timescale of 
that burst must have been of order 110 Myr to account 
for the observed numbers of such bursts: 
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Here, ni,urst represents the number density of objects 
with Alog SFR > 0.5, and Unew quiescent equals the 
number of newly quenched (A log(S'i^i?) < — 1) systems 
since z — 1.5 (which is most of them, given the rela- 
tively low number density at z > 1.5). The time interval 
probed in the z ^ 1 bin spans Aio.5<z<i.5 = 4.2 Gyr. 
For the z '-^ 2 bin, similar arguments lead to a value of 
^ 70 Myr. To determine these timescales, we restricted 
ourselves to the deepest fields, GOODS-South and UDS, 
to avoid any incompleteness effects. Our WFC3 iJieo- 
selected samples in these fields do not suffer from incom- 
pleteness above 10^" M©, even for maximally old, high 
M/ L systems in our highest redshift bin (see Appendix 
B.l). Increasing the number statistics by pushing down 
to 5 X 10^ Mq, where we are still highly complete, we 
obtain timescales of 170 Myr and 120 Myr at z ~ 1 and 
z ~ 2 respectively. The derived timescales should be con- 
sidered as crude upper limits. First, uncertainties in SFR 
estimates may lead to a contamination in the high-SFR 
tail by galaxies with intrinsically lower SFR. Second, in 
the case that not all excursions to high SFRs are followed 
by quenching, the inferred timescale of bursts associated 
with a quenching event will be shortened further. Over- 
all, the timescales obtained are not inconsistent with 
those expected from merger-induced star formation on 
the basis of hydrodynamical simulations (Mihos & Hern- 
quist 1996; Cox et al. 2008), although there may be room 
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for other evolutionary paths to reach the quiescent phase 
as well. 

Our observations of the structural MS and the qui- 
escent population located below in Figure 1 reveal two 
aspects that any galaxy formation and quenching model 
should account for. First, we find spheroid-like (n ^ 4) 
quiescent galaxies to be present over a wide range of 
stellar masses (1.5 - 2 dex) at all redshifts, down to at 
least 10"'^° Mq. Second, galaxies with a large spread in 
star formation activities co-exist at the same stellar mass 
over a mass range of an order of magnitude (roughly 
10 < log(M) < 11). If all galaxies at a given stellar 
mass are central galaxies and reside in an equally mas- 
sive dark matter halo, this would imply that quenching is 
not a simple step function of halo mass. Carrying out a 
detailed abundance matching analysis between galaxies 
at < z < 1 and their halo hosts, Conroy & Wechsler 

(2009) also do not find evidence for a sharp characteristic 
halo mass at which star formation truncates. If on the 
other hand halo quenching is responsible for the transi- 
tion from active to passive systems, this would imply that 
quiescent galaxies live in higher mass halos (i.e., have 
a smaller baryon fraction) than SFGs of similar stellar 
mass. 

Finally, a third quenching mechanism that can con- 
tribute to the build-up of the quiescent population is en- 
vironmental (satellite) quenching. Studying the role of 
environment and its impact on the morphological prop- 
erties of (part of) the quenched population, is beyond 
the scope of this paper. We do however note that Y. J. 
Peng et al. (2010) argue this physical process only be- 
comes relevant at late times (z < 0.5) and low masses 
(logM < 10.5). Based on an empirical model that suc- 
cessfully reproduces the stellar mass function of star- 
forming and passive galaxies in SDSS, Y. J. Peng et al. 

(2010) attribute the quenching at larger lookback times 
to a combination of a mass-dependent component (dom- 
inating at the high-mass end), and a mass- independent 
merger quenching component. Together, these processes 
are claimed to be responsible for the double Schechter 
mass function of passive galaxies. It is important to point 
out that their 'mass quenching' is probabilistic in nature, 
rather than a step function, with a rate proportional to 
the SFR, and an underlying physical origin that may 
well be associated with AGN or star formation feedback 
rather than a halo-scale phenomenon. 

5.2. A Physical Model for IR/UV Ratios 

In Sections 4.1 and 4.2, we presented observational re- 
sults on how the surface density and obscuration of star 
formation depends on the position of galaxies in the SFR- 
Mass diagram. Here, we investigate whether those two 
observations can be understood in the context of one self- 
consistent physical model that relates the amount of re- 
processed emission to the column of obscuring material, 
and indirectly to 'Esfr- The aim of such a model is to 
explain qualitatively the difference in shape of iso-'Esfr 
and iso-SFRin/SFRuY regions in the SFR-Mass plane 
(Figure 4 versus Figure 6), and to reproduce quantita- 
tively the observed SFRm/SFRuy ratio starting from 
"^SFR, with the help of a few fundamental, and where 
possible empirically calibrated, laws. 

Our model involves three steps: a conversion from SFR 
surface density to gas surface density, a conversion from 



gas column to dust column, and a geometry-dependent 
translation from dust column to IR-over-UV ratio (see 
top panel of Figure 7) . The first step boils down to an 
application of the Kennicutt-Schmidt (KS) star forma- 
tion law (Kennicutt 1998). Recently, Genzel et al. (2010) 
made the case for a KS law of the form 

\ogiT. sFRiM^yr-'^kpc-^]) = 1.171og(S„oj gas[MQPc-^]) 

-3.48 (4) 

on which both nearby and high-redshift main sequence 
SFGs lie. Major mergers follow a relation of similar slope 
according to these authors, but may be offset in nor- 
malization, so that a higher level of star formation is 
obtained for a given reservoir of molecular gas. For sim- 
plicity, we assume Equation 4 holds for all galaxies in our 
sample, at low and high redshift. We apply Equation 4 
to our observed "Esfr, account for a 36% mass correc- 
tion for Helium to obtain the H2 surface mass density, 
and add a fiducial plateau of HI surface mass density 
Y^Hi = 4.5 Mqpc'^ (Krumholz, McKee & Tumlinson 
2009) to obtain the total hydrogen column. In most 
cases, the contribution of HI to the total gas column is 
negligible. 

The second step of our model accounts for the fact 
that not all galaxies have the same metallicity. This is 
relevant, because the limited number of cases in which 
metallicities as well as dust-to-gas ratios have directly 
been observed (Milky Way, Small and Large Magellanic 
Cloud), illustrate that it is not the dust-to-gas ratio, 
but rather the dust-to-metal ratio that remains relatively 
constant between galaxies. We adopt a gas-to-dust ra- 
tio equal to that of the Milky Way {Nh/Av)mwg = 
1.87 X 10^1 cm-'^ mag-^ (Bohlin et al. 1978), with a 
linear scaling accounting for gas-phase metallicities devi- 
ating from solar: 



Nh/Av = {Z/Zq)-\Nh/Av)mwg- 



(5) 



Even in the era of multi-object near-infrared spectro- 
graphs, spectral metallicity tracers of high-redshift galax- 
ies, such as [Nil]/ Ha are notoriously difficult to ob- 
tain for large samples of galaxies, let alone for the sam- 
ple size considered here. We therefore adopt an empir- 
ically calibrated functional form to compute the metal- 
licity on the basis of known stellar population proper- 
ties. Rather than assuming an evolution of the mass- 
metallicity relation with redshift Z{M,z) (e.g., Erb et 
al. 2006; Maiolino et al. 2008; Mannucci et al. 2009; 
Buschkamp et al. in prep), we parametrize the metal- 
licity as a non-evolving function of position in the SFR- 
Mass diagram: Z{SFR,M). Using the SDSS spectral 
database, Mannucci et al. (2010) established such a fun- 
damental surface in SFR-Mass-Metallicity space. Fur- 
thermore, these authors find that galaxies out to z ^ 2.5 
lie on the same surface as do nearby systems. The ob- 
served evolution of the mass-mctallicity relation stems 
in this picture merely from the fact that galaxies with 
higher SFRs are more abundantly present at high red- 
shift. That higher SFR systems at a given mass have 
lower metallicities is interpreted as evidence for recent 
gas infall that simultaneously fuels the formation of new 
stars, and dilutes the metal-content of the interstellar 
medium. 
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Fig. 7. — A comparison between the observed SFRju/SFRjjy ratio and that inferred from the observed T,sfr using a simple physical 
model that takes into account the SF law, metallicity dependence, and geometric distribution of dust and stars. The solid and dashed 
contours mark the 25th and 75th percentiles of the galaxy distribution respectively. While at 2 ~ 0.1 the model agrees with the observations, 
towards higher redshift a relatively larger amount of UV radiation is observed than expected from the model. This may infer shortcomings 
related to the assumed geometry, or alternatively SF law and/or metallicity dependence. 



Even in the local universe, measurements of the gas- 
phase metallicity (and therefore the normalization of the 
mass-metallicity relationship) vary by factors of ^ 3 
depending on the tracer and calibration used (Kewley 
et al. 2008). Given this systematic uncertainty, we 
choose to rcnormalize the Mannucci et al. (2010) re- 
lation so that it reaches a plateau of solar metallicity 
{log{0/H) + l2 = 8.66, Asplund et al. 2004) at the high- 
mass end, so that for galaxies with Milky Way properties 
(~ 10" Mq, 1 MQ/yr), the Milky Way dust-to-gas ratio 
applies: 

\og{0/H) + 12 = 8.49 + 0.47(mo.32 - 10) if M0.32 < 10.2 
8.66 if /io.32 > 10.5 (6) 

where /io.32 = log(M) - 0.321og(5Fi?). For 10.2 < 
Mo. 32 < 10.5, the smooth polynomial form from Man- 
nucci et al. (2010, equation 4) is applied, downscalcd by 
the same factor. 

The scaling of the column with Z{SFR,M) accounts 
at least in part for the difference in shape of the Yisfr 
and SFRir/SFRuv diagrams (Figure 4 and 6 respec- 
tively). While along the MS little evolution in Y,sfr is 
observed, gas-phase metallicities are expected to increase 
as we march up the MS. This leads to an increase in the 
dust column, and consequently of the reprocessed frac- 
tion of the star formation towards the high-mass end of 
the MS, as observed in Figure 6. In other words, our 
observations do not probe metallicities directly, but are 
consistent with the presence of a mass-metallicity and/or 
mass-SFR- metallicity relation. ^^ 

^® The assumption of solar metallicity in determining the stellar 
masses (see Section 2.2.2) is in principle inconsistent with the gas- 
phase metallicity scalings applied in our model. As a sanity check, 
we verified that all conclusions drawn in this paper remain valid if 
we were to adopt a stellar metallicity equal to the gas-phase metal- 
licity according to Mannucci et al. (2010); Z* = Z{SFR,M{Zt)). 



Using the Calzetti et al. (2000) attenuation law in com- 
bination with Equation 5, we now compute the optical 
depth r28oo at the rest-frame 2800A wavelength where 
we defined SFRjjv' 



T26m = (1.79Av)/1.086 



(7) 



The third and final step of our model relates the opti- 
cal depth to the ratio of reprocessed over unobscured star 
formation. For the geometry of a uniform foreground 
screen, i.e., where all stars are hiding behind the total 
dust column, the intrinsic UV emission from young stars 
is reduced by a factor e""^. We quickly dismissed such 
a geometry, as it is physically implausible, and leads to 
inferred SFRjji/ SFRnv that exceed the observed val- 
ues by several orders of magnitude. Clearly, UV emission 
from young stars is observed, and assuming a more real- 
istic mixture of dust and stars in our model can account 
for this. We choose to adopt the simplest possible model, 
in which the emitting sources and the dust are homo- 
geneously mixed. In such a configuration, the intrinsic 
emission is attenuated by a factor (1 — e~'^)/r (McLeod 
et al. 1993; Forster Schreiber et al. 2001). Rewritten to 
a functional form for SFRju/SFRuv, this becomes: 



SFRijii SFR, 



^2800 



1 



UV 



1 



-''"2800 



(8) 



We present the comparison of our model results to the 
observed values of SFRm/ SFRjjv in Figure 7. Solid 
and dashed contours enclose 25% and 75% of the galaxy 
distribution respectively. In the z '-^ 0.1 bin, passive 
galaxies, for which our SFRin/SFRjjv is ill-defined (see 

The effect of such a scaling is to slightly skew the SFR-Mass dia- 
gram, such that generally lower mass galaxies (and at a given mass 
particularly those with higher SFRs) shift to masses that are lower 
by to -0.2 dex. At metallicities of a half solar, the stellar masses 
decrease by a few percent in the median, while the median decrease 
amounts to -0.1 dex at a quarter of solar. 
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Section 4.2) were excluded from the analysis. For nearby 
galaxies, and even at z ^ 1, we find a surprisingly good 
agreement between model and observations. A slight 
bias is present at z ~ 1, in the sense that the inferred 
SFRju/SFRifv for the bulk of the galaxies exceeds the 
observed values by a factor of ^ 2. This systematic off- 
set increases further toward higher redshift, amounting 
to an order of magnitude at z ^ 2. 

We speculate briefly on the physical origin of this 
trend. A refined model should account for an increased 
amount of UV emission escaping unhindered from high- 
redshift galaxies. This situation is reminiscent of Adel- 
berger & Steidel (2000), who found that for galaxies of 
a given SFR, at high redshift a larger proportion of that 
star formation is emitted in the UV. Daddi et al. (2007) 
also discussed that z ~ 2 galaxies are more " transparent" 
to UV radiation than local galaxies with similar levels of 
star-forming activity. The offsets found in the z ~ 1 and 
z ^ 2 panels of Figure 7 could signal shortcomings in 
any of the three steps of our simple model. Since the 
first two steps (star formation law, and metallicity scal- 
ing) are empirically motivated, and the third (geometry) 
is not, it makes sense to address the latter first. 

Given that star formation, dust production, and the 
physical processes (such as winds) that determine the 
distribution of the ISM are intimately connected, the 
idealized model of a homogeneous mixture of dust and 
(young) stars will likely break down at some level. If 
high-redshift galaxies are more patchy than their local 
counterparts, this would naturally explain why a rela- 
tively larger fraction of the intrinsic emission can escape 
the galaxies unattenuated. In contrast to galaxies today, 
the UV and IR emission of high-redshift systems would 
in that case more often be associated with physically dis- 
tinct regions within a galaxy. If such a patchy configu- 
ration would indeed be more representative for z ^ 2 
galaxies, we would expect this to be revealed in terms of 
rich internal color variations in multi-wavelength high- 
resolution imaging. Such internal color variations, and 
their interpretation in terms of spatially resolved stel- 
lar populations will be the focus of a forthcoming pa- 
per (Wuyts et al. in prep). First steps in this direction 
have been taken through visual classification (Cameron 
et al. 2010) and quantitative internal color dispersion 
methods (Bond et al. 2011) to compare ACS and WFC3 
morphologies of distant galaxies in the ERS and HUDF 
regions. 

Alternatively, the results of our model comparison 
may signal deviations from the adopted star formation 
law (Equation 4). A tuning of the KS law to en- 
force an agreement between the modeled and observed 
SFRjji/ SFRijv would require a higher normalization 
and/or a steeper slope than 1.17, so that a given Tjsfr 
corresponds to a smaller column of obscuring material. 
Genzel et al. (2010) demonstrate that for a subsample 
of galaxies (those identified as mergers) a higher normal- 
ization is indeed appropriate. Depending on their con- 
tribution to the overall galaxy population, this may lead 
to a somewhat steeper slope of the effective KS relation 
fitted to the combined sample of normal SFGs and merg- 
ers. The change in slope required to bring our model pre- 
diction in agreement with the observations (keeping the 
zeropoint fixed) is substantial: a slope of 1.35 at z ~ 1 
and 1.75 at z ^ 2. KS slopes of this magnitude have been 



reported in the literature (e.g., Bouche et al. 2007), but 
in combination with a lower zeropoint of the relation. If 
instead we keep the slope fixed to 1.17, a zeropoint in- 
crease to -3.2 at z ~ 1 and -2.5 at z ~ 2 would be required 
to obtain agreement between model and observations. A 
large fraction of z ~ 2 galaxies for which gas surface 
densities have been measured do lie on a relation with 
slope 1.17 and zeropoint -2.44, but that is because due 
to sensitivity limitations these samples have often been 
biased to high-SFR outliers. Recent measurements for 
z ^ 2 galaxies on the main sequence (i.e., representative 
for the bulk of the z ~ 2 population) are consistent with 
the local relation for normal SFGs (Genzel et al. 2010; 
Daddi et al. 2010), and fitted by themselves give a zero- 
point of -3.11 (and slope of 1.17). The required change 
therefore seems to be in tension with the observational 
constraints at z ^ 2 available to date. 

Finally, a lower dust column for a given Yisfr could 
also be obtained if our metallicities of high-redshift galax- 
ies are overestimated. If this is the case, this would imply 
that gas-phase metallicity is determined by three param- 
eters {Z(SFR, M, z)) rather than a non-evolving depen- 
dence on SFR and stellar mass. We stress that current 
observations out to z ^ 2.5 do not require such an evolu- 
tion with redshift of the SFR-Mass-Metallicity plane, but 
the high-redshift samples presently show a large scatter 
(Mannucci et al. 2010). Moreover, while we are still fac- 
ing large systematic uncertainties in metallicity measure- 
ments of even local galaxies, we already took the liberty 
to rescale the Mannucci et al. (2010) relation to reach a 
plateau of solar metallicity at the high-mass end, which 
is among the lower end of calibrations presented in the 
literature (Kewley et al. 2008). Vice versa, if we were to 
adopt the Mannucci et al. (2010) relation without rescal- 
ing (i.e., reaching a plateau of 12 -I- \og{0/H) = 9.07, or 
2.6 times solar), the model prediction of SFRjji/SFRuv 
would exceed the observed values by a factor 5, 10, and 
60 at z ~ 0.1, z ~ 1, and z ~ 2. 

We conclude that a simple model, combining a mixed 
dust-star geometry with an empirical star formation law 
and metallicity scaling, is successful in describing the 
obscuration and surface density of star formation in a 
self-consistent picture. Towards higher redshifts, an ad- 
justment of this model is desired, allowing more of the 
UV emission to escape. This may be due to the fact that 
galaxies at higher lookback times feature more patchy 
geometries. 

6. SUMMARY 

We use statistically significant samples of unprece- 
dented size (~ 640000 galaxies at z - 0.1, - 130000 at 
z '^ 1, and ^ 36000 at z ^ 2) to analyze how the struc- 
ture and mode of star formation of galaxies at different 
lookback times depends on their position in the SFR- 
Mass diagram. To this end, we combine 4 deep lookback 
surveys (COSMOS, UDS, GOODS-N, GOODS-S) with 
high-resolution imaging and rich multi-wavelength data 
sets reaching from the optical and near-infrared to mid- 
and far-infrared wavelengths. Value-added SDSS cata- 
logs are complemented with GALEX data to investigate 
equivalent galaxy properties in the nearby universe. 

We find remarkable similarities over 10 Gyr of lookback 
time. E.g., the relation between surface brightness profile 
shape on the one hand and the position of a galaxy with 
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respect to the main sequence of star formation on the 
other hand is strikingly similar at all observed epochs. 
Phrased differently, we find evidence for a correlation 
between galaxy structure and stellar populations (i.e., 
a Hubble sequence) to be in place already three billion 
years after the Big Bang. The fact that a structurally 
distinct quiescent population below the main sequence is 
already present since z ~ 2.5 implies that the quench- 
ing mechanism responsible for the shutdown of star for- 
mation must be universally present already since that 
epoch, and must be causally connected to the morpho- 
logical transition as well (or vice versa, see Martig et al. 
2009 who analyze a morphological quenching scenario). 

Low- and high-rcdshift galaxies also seem to follow a 
qualitatively similar behavior in terms of star formation 
surface density and its dependence on the specific SFR 
(see also Elbaz et al. 2011), and the variation in ob- 
scuration (parametrized by SFRm/SFRuy) along and 
across the MS. They differ in the sense that the high 
end of the relation between ^sfr and SFR/M is more 
abundantly populated at larger lookback times. 

Another aspect in which low- and high-redshift galax- 
ies differ is their typical size. Measured at fixed 
[SFR, M] , galaxies on the MS today experienced a slower 
growth over time (~ (1-1- 2;)^'^'^) than the quiescent pop- 
ulation (~ (l-|-z)~^-^). These growth rates are consistent 
with earlier work based on ground-based imaging (Tru- 
jillo et al. 2006; Franx et al. 2008; Toft et al. 2009) and 
recent WFC3 measurements by Nagy et al. (2011).^^ A 
simple mass-size relation is inadequate to describe the 
size variations of the entire galaxy population. At all 
epochs, star-forming galaxies on the MS are the largest 
among galaxies of a given mass. 

The profile shape of galaxies, quantified by the Sersic 
index n, is not a monotonic function of the level of star 
formation. Instead, at a given stellar mass a reversal to- 
wards higher n is observed in the high-SFR tail of the 
distribution. Both our low- and high-redshift samples 
show evidence of such a relation. This trend may hint at 
an evolutionary connection between the high-SFR out- 
liers, which are in the process of building up a nuclear 
concentration of mass, and the quiescent galaxy popula- 
tion, which is characterized by cuspy surface brightness 
profiles, as proposed by the merger paradigm. 

The variation of Ssf_r across the SFR-Mass diagram 
shows a different behavior than that of SFRm/ SFRuv 
While iso-Tjsfr regions run diagonal, following approx- 
imately lines of constant specific SFR (see also Schimi- 
novich et al. 2007 for z - 0.1), the SFRir/SFRuv 
ratio increases as we move up along the MS (i.e., at fixed 
specific SFR) . This strongly suggests variations in metal- 
licity along the main sequence, as expected from stud- 



ies of the mass-metallicity (Erb et al. 2006; Buschkamp 
et al. in prep) and SFR-mass-metallicity (Mannucci et 
al. 2010) relation. An increase in gas-phase metallic- 
ity towards the high-mass end of the MS would cause 
the obscuring column to increase, even if the gas column 
remains the same. 

Using a simple model involving the observed 
Kennicutt-Schmidt star formation law (Genzel et al. 
2010) and an empirically calibrated metallicity scaling 
(Mannucci et al. 2010) in combination with an as- 
sumed dust-star geometry, we compare the observed 
SFRir/SFRjjv ratios to those inferred from the ob- 
served YisFR- Our model is successful at qualitatively 
relating the two properties, and explaining why their iso- 
contours in the SFR-Mass plane are of different shape. 
For a metallicity scaling that flattens at solar values at 
the high-mass end, we find a good quantitative corre- 
spondence at 2; ~ 0.1 and 2^1 when assuming a homo- 
geneous mixture of dust and stars. At higher redshifts, 
progressively more patchy dust-star geometries appear 
to be required in order to have sufficient UV emission 
escaping to match the observed SFRm/SFRiJv 
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^'' Nagy et al. (2011) find similar re{z)/re{z = 0) as star-forming 
galaxies in our sample that satisfy the same selection criteria. Their 



steeper growth rate (r^ (l + z) 



within 1.5 < 2 < 3) stems from 



the fact that they do not anchor the size evolution to z = 0, i.e., 
they effectively use a different parametrization. 



APPENDIX A GALAXY PROPERTIES ACROSS THE MAIN SEQUENCE 



Similar to Figure 2 in Section 3.1, we here present the dependence of the effective radius r^, surface density of star 
formation 'Esfr, and obscuration SFRiji/ SFRjjv of galaxies with 9.7 < logM < 11 as a function of deviation from 
the MS (Figure 8). The scatter at a given /S.\og{SFR) {grey polygons) stems from a combination of uncertainties in 
the derived parameters and intrinsic variations within the galaxy population of a given specific SFR. Especially the 
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relations reflecting the mode of star formation (Ssf_r and SFRm/SFRuv) with respect to the MS are remarkably 
tight. 
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Idem as Figure 2, but for log(re), log(Sg^jj), and log{S F Rj n/ S F Ruv) instead of Sersic index. 



APPENDIX B UNCERTAINTIES AND COMPLETENESS 



We follow two approaches to address the robustness of the trends in SFR-Mass space described in this paper. First, 
we compare empirically the results for individual fields (B.l) and for data sets of different depth and wavelength 
in a given field (B.2). Second, we use Monte Carlo simulations (B.3) to propagate different sources of uncertainty 
(photometry, redshifts, SFRs, structural parameters, ...) and assess their impact on the median properties of the 
overall galaxy population. While specific measurements on individual objects may sometimes be subject to significant 
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Table 1. Overview Deep Lookback Surveys 



Field 



Area 
(dcg2) 



Filter, 



morph 



Image Depth* 
(AB mag, Scr) 



Sample Depth'' 
(AB mag) 



No.5 



Nl.5 



COSMOS 


1.480 


-^814 


27.2 


25.0 


106080 


21430 


UDS 


0.056 


^160 


26.7 


26.7 


10443 


6796 


GOODS-S 


0.041 


^160 


27.0 


27.0 


7008 


3973 


GOODS-N 


0.042 


2850 


27.6 


26.8 


8797 


3450 



''Point source depth of the image on which the morphological analysis was performed. 

''Magnitude (in i, Hiqq, Hiqq and zgso for COSMOS, UDS, GOODS-S and GOODS-N respectively) down to which galaxies were included 
in our sample. 
'^Sample size in the 0.5 < 2 < 1.5 and 1.5 < 2 < 2.5 redshift intervals. 



uncertainties, the conclusions drawn in this paper are much less affected as they are based on median properties over 
large numbers of galaxies, and the level of systematic biases is small. 

B.l Field-to-Field Variations and Completeness 

Our results at intermediate and high redshift are based on 4 lookback surveys, each with a wealth of multi- wavelength 
data. They complement each other in area and depth, and differ in some cases in the selection band of the multi- 
wavelength catalog as well as the resolution and waveband of available HST imaging (see Table 1). A comparison of the 
results obtained for each of the fields individually therefore serves as an empirical robustness check that encompasses 
field-to-field variations, the impact of imaging depth on morphological measurements and completeness within a given 
[SFR, M] bin, resolution, and morphological k-corrections, if present. In Figure 9, we demonstrate that the median 
Joehavior of n, re, T,sfr and SFRm/SFRuv looks strikingly similar for our samples in COSMOS, UDS, GOODS-S and 
GOODS-N, despite inhomogeneities between the data sets. In cases where the patterns look somewhat less smooth, 
this is to be expected from the smaller number statistics. From the similarity of the T,sfr and SFRiji/ SFRjjv 
diagrams, it follows that also our simple model and interpretation in terms of patchier dust geometries at high redshift 
holds for each of the fields separately. 

The region of SFR-Mass space covered by each of the fields in Figure 9, reflects differences in depth. We expand on 
this in Figure 10, where the color-coding illustrates the estimated completeness of our sample. For the COSMOS and 
GOODS-N fields, the completeness in each [SFR, M] bin was defined empirically as the fraction of UDS and GOODS-S 
galaxies that would enter the i- and z-band limited samples for those fields respectively. For the deep WFC3 _ff-band 
selected samples in UDS and GOODS-S, we estimated the completeness indirectly. For a given [SFR, M] bin, we 
dimmed the SEDs of galaxies with similar specific SFR, but higher SFR and mass (i.e., located in a highly complete 
region of the diagram) . As completeness estimate, we adopted the fraction of galaxies downscaled to the [SFR, M] 
bin of interest that would still be sufficiently bright to enter the sample. The latter is best interpreted as an upper 
bound on the completeness. As extinction is known to scale with the star formation activity (see, e.g., the increased 
obscuration as one moves along the MS toward the high-mass end in Figure 6), our simple approach of downscaling 
the SED by the same factor as the SFR and mass, without altering the SED shape, likely leads to a completeness 
estimate that is too conservative. 

Figure 10 illustrates that the completeness level of our COSMOS sample is lower than that of the samples extracted 
from the other fields. As is evident from Figure 9, COSMOS therefore gets zero weight at the lowest SFRs and masses 
probed in the z ~ 2 bin. In [SFR, M] bins that do contain COSMOS galaxies, two tests indicate that the somewhat 
lower completeness level does not lead to significant biases in the median morphology or mode of star formation. 
First, the trends in COSMOS do not differ from those observed in the other fields (Figure 9). Second, we applied the 
same i-band limit adopted for COSMOS to the highly complete WFC3-selected catalog in UDS. We verified that this 
does not change the median morphological and mode of star formation parameters in [SFR, M] bins that still contain 
sources after imposing the i-band cut. 

B.2 Depth and Wavelength Dependence of Morphological Measurements 

While boosting confidence in the robustness of our results, the field-to-field comparison in B.l did not discriminate 
between different aspects of inhomogeneity between the data sets. Using the GOODS-S field as a test case, we now 
isolate the effect of depth, and the effect of waveband on which the surface brightness profile fitting is performed. 
The top panels of Figure 11 illustrate how the morphological parameters change when measured on the presently 
available full depth mosaic of the CANDELS-Deep region, and on a version of the mosaic built with only the first half 
of the exposures (0.35 mag shallower, similar in depth to the CANDELS-Wide data used for UDS). Unsurprisingly, 
the scatter in Alog(re) is observed to increase toward fainter magnitudes. However, down to the faintest magnitudes 
we find no evidence for a systematic bias in the size measurement. This implies that population studies which describe 
the median size, such as presented in this paper, can safely probe to fainter levels than would be considered reliable on 
an individual object level. Haussler et al. (2007) arrive at the same qualitative conclusion when inserting simulated 
n ~ 1 and n — 4 galaxies in mock ACS images, and recovering the morphological parameters using GALFIT. For 
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Fig. 9. — Prom top to bottom: SFR versus mass diagrams of the 4 deep fields separately, color-coded by median Sersic index, size, surface 
density of star formation, and SFRjji/ SFRjjy ratio. Despite differences in depth, and waveband used to perform the morphological 
analysis, the patterns in SFR-mass space are consistent in all fields, boosting confidence in the robustness of our results. 



similar tests aimed specifically at the CANDELS data sets, we defer the reader to van der Wei et al. (in prep). It is 
relevant to note that our conclusion may not necessarily hold for other methods to determine r,. than the GALFIT 
code used in our analysis (see also Haussler et al. 2007; Sargent et al. 2007). 

The difference An between the Sersic index measured on the 2- and 4-epoch mosaics also exhibits an increased 
scatter toward fainter magnitudes, but again the median An remains consistent with zero. In detail, however, a trend 
is noted such that shallow profile shapes are on average recovered correctly down to the faintest magnitudes, while 
the Sersic index of intrinsically cuspy systems shows a systematic deviation at the faint end.^® In relative terms, the 
underestimate for faint n > A galaxies corresponds to a median An/n « —0.4. In SFR-Mass space, these galaxies 
are predominantly located at log(M) < 9 and log(S'_Fi?) < 0.5. If they were the dominant population in that area 
of SFR-Mass space, they could potentially bias the median n to low values. Our deepest, 4 epoch GOODS-S data 
suggests however that they only account for less than 10% of the galaxy population in the low mass and low SFR 
regime. 

The bottom panels of Figure 11 compare measurements of n and re carried out on the WFC3 i?i60 mosaic, and on 
the ACS ^850 mosaic. We plot the deviation (negative where the zs5o measurement yields a smaller value) as a function 
of redshift to facilitate the interpretation in terms of morphological k-corrections in the rest-frame. Variations in the 
extent of the zsso- and -ffigo-band surface brightness profiles are small compared to the amplitude of the size variation 
across the SFR-Mass diagram. Median offsets are limited to the few percent level for r^, and of order ~ 10% for the 



^^ We note that deviations An of order -2 or lower can by definition not be reached for disk-dominated systems because of the allowed 
range 0.2 < n < 8. 
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Fig. 9. — continued 
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Fig. 10.— Completeness of our high-redshift samples in COSMOS, UDS, GOODS-South, and GOODS-North as a function of SFR and 
mass. 
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Fig. 11. — Comparison between different GALFIT runs in GOODS-South. Top panels: Deviation in haff-light radius and Sersic index 
measured on WFC3 Hieo mosaics of 2 epochs versus 4 epochs depth as a function of i?i60 magnitude. Bottom panels: Deviation in 
half-light radius and Sersic index measured on ACS zg^g versus WFC3 J^ieo imaging as a function of redshift. Horizontal red lines indicate 
the running median, while vertical red lines mark the central 50th percentile. Systematic deviations from zero are mostly limited to a few 
percent, much smaller than the variation of morphological parameters across the SFR-Mass plane. Offsets in Sersic index at the 10% level 
are present between zg^o and Higg, such that the surface brightness profiles are slightly cuspier in the longer wavelength band. 



Sersic index (such that the longer wavelength surface brightness profiles tend to be cuspier). OveraU, the systematic 
deviations in n are small compared to the amplitude of trends with profile shape across the SFR-Mass plane, although 
we caution that the random scatter in n due to morphological k-corrections becomes substantial for systems with 
n > 4 in the i^ieo band. 

B.3 Propagating Uncertainties with Monte Carlo Simulations 

We ran extensive Monte Carlo (MC) simulations to estimate the effects of photometric and photometric redshift 
uncertainties on our analysis. First, we selected a representative subset of galaxies in each of the four deep fields. 
We sample four redshift intervals (0.5 < 2; < 1.0, 1.0 < 2; < 1.5, 1.5 < z < 2.0, and 2.0 < z < 2.5), drawing for 
each one three galaxies from each [SFR^ M] bin, spanning the lower, middle and upper fertile in signal-to-noise. For 
the resulting subsample of 7480 galaxies, we then compiled one hundred mock photometric catalogs. The flux-points 
in each optical-to-IR band were drawn randomly from a Gaussian with the measured flux as the mean and its error 
as the standard deviation. Next, we subject each mock catalog to the same fitting procedures to derive photometric 
redshifts, stellar masses, and SFRs, and also propagate the newly obtained redshifts into the conversion from angular 
to physical half-light radii. Together, the 100 mock realizations provide us with an uncertainty distribution on (and 
often correlated between) each of the parameters that enter our analysis. 

In order to apply the MC results, we associate each galaxy from our full sample with a galaxy in the MC subsample 
that lies in the same field, and occupies the same redshift, SFR, mass and signal-to-noise bin. The deviations with 
respect to the actual measurements of redshift, mass, SFR, and r^ are drawn from one of the 100 mock realizations of 
the associated MC target, and applied to the galaxy of interest. The resulting mock population is then again divided 
in [SFR, M] bins and we compute the median of the property of interest (n, r^, T^sfr, SFRjr/ SFRuv), exactly as 
we did for our analysis of the real measurements. In principle, the deviations in derived parameters (which unlike the 
photometry are not assumed or forced to be symmetric) should be applied to error-free measurements. Since the latter 
are not available, our test serves as a conservative estimate of how important photometric errors are in our analysis. 

Upon visual inspection of the MC-ed equivalents of Figures 1 to 7, we find all the trends described in this paper 
to be robust against photometric uncertainties. However, the Zphot uncertainty distribution as inferred from running 
EAZY on the 100 perturbed photometric catalogs, is likely too narrow. We find evidence for this by running a similar 
set of 100 simulations on the full spectroscopically confirmed sample. The spectroscopic redshift lies within the formal 
Icr confidence interval derived from the MC simulation in only about 30% of the cases, much less than the expectation 
value of 68% were photometric errors the only source of uncertainty in deriving the redshift. 

In a second iteration, we therefore base deviations in Zphot on the empirical comparison between photometric redshifts 
and the available spectroscopic redshifts. Photometric and photometric redshift uncertainties are then propagated to 
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Fig. 12. — Monte Carlo simulated SFR vs mass diagrams at ^ ~ 1 and 2 ~ 2, constructed by perturbing the measured photometry and 
photometric redshifts according to their uncertainties and propagating to the other parameters of interest. The characteristic patterns in 
SFR-Mass space remain robustly present in the MC simulation. 

the other parameters as before. Here, it is important to account for the fact that the spectroscopic sample is not 
randomly drawn from the overall galaxy sample, but instead biased towards bright magnitudes. Following Scarlata 
et al. (2007), we quantified the photometric redshift uncertainty at a given magnitude semi-empirically by dimming 
the SEDs of spectroscopically confirmed galaxies to that magnitude (perturbing the photometry accordingly), and 
deriving the resulting Az/(1 + z) distribution. As expected, this yields a broader distribution, at the faintest levels 
by up to a factor of 5 compared to the normalized median absolute deviation of the undimmed spectroscopic sample. 
However, even at the faintest levels, no significant systematic offsets in Zphot are revealed. 

Figure 12 presents the MC-ed equivalents of Figures 1 to 6, where photometric errors and semi-empirical photometric 
redshift uncertainties have been applied and propagated in the analysis. We considered the distribution in Ap = pMC — 
ptruej where p stands for the median structural (n, re) or mode of star formation (Esfr, SFRjh/ SFRuv) property 
in a given [SFR,M] cell. Systematic deviations from zero are limited to < 0.03 dex. At 0.5 < z < 1.5, the scatter 
of the Ap distribution is modest, amounting to 0.04, 0.04, 0.07, and 0.12 dex in n, r^, '^sfb. and SFRibJ SFRijv 
respectively. At 1.5 < z < 2.5, the scatter in Ap doubles. We conclude that, when following this more conservative 
approach, we still find that the patterns in SFR-Mass space described in this paper are robust against photometric 
and photometric redshift uncertainties. 

In a third and final iteration, we now also include random uncertainties in mass, SFR, r^., and n that do not stem 
from photometric or photometric redshift uncertainties. To this end, we adopt fiducial errors of 0.2 dex in each of 
these parameters. This crudely reflects our ability to recover the intrinsic surface brightness profile shape (B.2), as well 
as imperfections with which observables are translated to physical parameters (see, e.g., Wuyts et al. 2009, 2011a). 
The scatter in Ap is boosted by ~ 50% compared to the MC simulation that only accounted for photometric and 
photometric redshift uncertainties. Nevertheless, the overall conclusions drawn in this paper remain valid. 
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